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POTENTIAL FLOW COMPUTER PROGRAM* 
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SUMMARY 


A potential flow computer program has been in use at the Arnold Engineering De- 
velopment Center (AEDC) for several years. This program has been used primarily as 
a tool for flow-field analysis in support of test activities in the transonic wind tunnels 
of the Propulsion Wind Tunnel Facility (PWT). Analyses have been made over a Mach 
number range from 0 to 0. 9 for a variety of configurations from aircraft to wind tun- 
nels, with excellent agreement between calculated flow fields and measured wind tunnel 
data. Analytical and experimental data for seven different flow analysis problems are 
presented 'n this paper. 
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INTRODUCTION 


The AEDC Three-Dimensional Potential Flow Computer Program (PFP) in the 
existing form was developed primarily as a result of the need to make calculations of 
the flow field in the vicinity of aircraft fuselages (typically at locations where aircraft 
inlets might be located). 'This need arose because of the support the theoretical flow- 
field calculations could lend to a research program carried out at AEDC to simulate the 
inlet flow fields in a wind-tunnel test of full/scale inlet/engine systems (refs. 1 and 2). 
Much of the computing capability that the PFP presently has resulted from these flow- 
field calculations which have as their primary variables the flow angularity (upwash and 
sidewash) over a y-z plane. After the initial solution of the velocity field for a given 
model attitude and Mach number is obtained, the upwash and sidewash can be deter- 
mined for any given point or over any grid desired. A new solution is required for each 
model attitude or Mach number. In addition to computing the upwash and sidewash, the 
PFP also computes the local Mach number, Cp, and flow streamlines. A computer 
plotting program has been written to supplement the PFP, and computer plots can be 
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obtained for most of the above parameters. A streamline can be traced from any point 
in the stream either upstream or downstream (or both). The capability to make a plot 
of the mathematical representation of the model geometry before running the complete 
program allows corrections to be made, if needed, with only a slight loss of computer 
time. A two-volume report (ref. 3) that details the program, modeling techniques, ap- 
plication, and verification has been published. 

This paper briefly outlines the PFP application to seven flow analysis problems in 
support of the transonic wind tunnels in the PWT at AEDC. 


SYMBOLS 

Values are given in both SI and U. S. Customary Units. The measurements and 
calculations were made in U. S. Customary Units. 

A angle of attack 

0^ lift coefficient 

C p pressure co efficient 

M m free- stream Mach number 

p/p 0 ratio of surface static pressure to free stream total pressure 
X coordinate along tunnel axis, positive downstream 

Y horizontal coordinate, sign as indicated 

Z vertical coordinate, sign as indicated 

a model angle of attack, deg, positive up 

0 model angle of yaw, deg, sign as indicated 

£ upwash, deg, positive up 

a sidewash, deg, positive as indicated 


PFP APPLICATION AND UTILIZATION 


The PFP at AEDC has been used primarily as a >ool for analysis of the flow in the 
far field. (Far field refers to a distance away from the analysis model surface equal 
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to, or greater than, the vortex spacing in the direction of flo\ \ ) The modeling tech- 
niques required for this type of analysis are presented in vc iume II of reference 3. 
However, work is underway to develop the modeling technique to allow accurate analysis 
of the surface pressure. Results of some of this continuing effort are reported in 
reference 4. All the flow problems presented here are of the far-field type. 


Flow Field Between Two Hollow Circular Cylinders 

The analysis of these cylinders was part of a research program in which the objec- 
tive was to create flow fields by some auxiliary method to simulate the flow entering a 
full-scale inlet/engine at high angles of attack and yaw. The device was to deflect (or 
induce) the flow upward as it passed between the inclined cylinders. A mathematical 
model of the cylinders is shown in figure 1. The last circumferential ray on each cyl- 
inder had trailing vortices that were trailed at an angle equal to one -half the cylinder 
pitch angle. A comparison between the theoretical and experimental flow angularity 
data is shown in figure 2. The theoretical results are shown as lines of constant flow 
angle (both upwash and sidewash), and the solid symbols show the relative location of 
experimental data with the magnitude of the measured angles indicated. The Mach num- 
ber at which these data were taken was 0. 9. It can be seen that the PFP overestimated 
the flow inclination angles by approximately 1°. 


Flow Field Around an Aircraft Fuselage 

Primary purpose of this analysis was to verify the results obtained from the PFP. 
Experimental flow-field data used for comparison with theory were available from wind- 
tunnel measurements made during the Tailor Mate test series. The objective of the 
wind-tunnel test was to determine the flow' field (upwash and sidewash) at a typical en- 
gine inlet fuselage location. The mathematical model of the fuselage configuration is 
shown in figure 3. The comparison between the predictions from the PFP and the wind- 
tunnel data for a pitch angle of 25° and a free-stream Mach number of 0. 9 are shown in 
figure 4. Here again, excellent agreement was obtained. 


Flow Field Under a Fuselage-Wing Configuration 

The purpose of this analysis was also for program verification; again experimental 
data obtained during the Tailor Mate studies were used. The fuselage-wing configura- 
tion was analyzed to compute the flow field under the wing at the wing-fuselage junction. 
The computer math model used in the analysis is shown in figure 5. A comparison be- 
tween the upwash and sidewash predictions and the experimental data for a Mach num- 
ber of 0. 9 and an angle of attack of 10° is shown in figure 6. Analytical and experi- 
mental data trends show excellent agreement, although the predicted data show some- 
what higher flow angularity gradients across the survey area than the measured data. 
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Streamtube Entering Inlet Behind Wing 


The purpose of this analysis was to determine the origin nf the streamtube entering 
the inlet in support of an inlet hot gas ingestion investigation. The comparison between 
experimental and theoretical data are shown in figure 7. These data were taken during 
a store separation study in an effort to verify the accuracy of the PFP to predict the 
correct flow field above and behind the wing. Data were taken at a Mach number of 0. 3 
with an angle of attack of 8°, and show excellent agreement between the experimental 
and theoretical values. The mathematical model and the predicted streamtube are 
shown in figure 8. The streamtube was determined by tracing streamlines from four 
locations beginning just upstream of the inlet and extending forw' rd to just upstream of 
the aircraft nose. A Mach number of 0. 3 and an angle of attack of 8° were also used 
for the streamtube analysis. This mathematical model is the largest analyzed to date, 
with 1559 loop vortices and 20 horseshoe vortices, and required approximately 4 hours 
run time on the AEDC IBM 370/165 computer. 


Inlet/Engine in Crosswind 

This analysis was made in support of a crosswind experiment conducted during an 
inlet study in the AEDC 16— ft (4. 88-m) Transonic Wind Tunnel (PWT-16T). The objec- 
tive of the analysis was to determine if a 0. 91-m-diameter (3-ft-diameter) crosswind 
simulator would adequately simulate the crosswind when used in conjunction with the in- 
let model, and to determine the position for the simulator to give best results. The 
theoretical analysis was made with the inlet/engine in an infinite crosswind. The math- 
ematical model included only a portion of the experimental model as shown in figure 9. 

A computer plot of the mathematical model is shown in figure 10. The engine ducts 
were closed on the downstream end and a negative source was located near the rear 
center of each engine duct to produce the correct inlet mass flow when that particular 
engine was in operation. Streamlines were traced upstream from near the four corners 
of the inlet, for each engine in operation, to determine the flow pattern of the air- 
stream entering the inlet. By tracing the streamlines, a fan position was deter- 
mined that would influence the inlet flow for all engine power settings and cross- 
wind velocities required. A typical flow pattern for the analysis is shown in figure 
11 for a crosswind velocity of 20.57 m/sec (67.5 ft/sec) with both engines operat- 
ing. 


Pressure Distribution in PWT-16T Contraction Section 

The objective of this analysis was to determine the pressure distribution 
along the bottom and side walls of the PWT-16T contraction section. Pressure dis- 
tributions were needed for use in a theoretical boundary-layer analysis of the wind- 
tunnel nozzle to support a test-section flow angularity study. Mathematical modeling 
used in the analysis is shown in figure 12. The flow in the test section area was speci- 
fied to give the pressure ratio desired for Mach number 0. 6. The analysis provided 



streamline information at a distance of 0. 305 m (1 ft) from the walls, and the calculated 
Cp was converted to p/p Q . Following the calculation of the theoretical pressure dis- 
tribution, the pressure distribution was experimentally measured in the contraction 
section. A comparison between the theoretical and experimental pressure distribution 
is shown in figure 13, with excellent agreement indicated. 


Strut Effects Analysis 

The objective of this analysis was to determine the strut effect corrections to mea- 
sured force and moment data for a slender winged vehicle with a mid-strut mount. The 
vehicle wing was located just forward of the strut. For this analysis the upwash angle 
was determined with the PFP for the body alone (fig. 14) and the body with strut (fig. 15). 
An incremental upwash angle was then determined at the wing location from these two 
sets of data. In this case the incremental values were negative because of the down flow 
around the strut. The incremental values along the wing location were averaged and the 
ACl correction calculated from the average angle -of-attack change. A comparison of 
the calculated corrections and those measured with a subscale model are shown in fig- 
ure 16. Excellent agreement is shown in both the trend with Mach number and the ab- 
solute values. 


CONCLUDING REMARKS 


The AEDC Potential Flow Program is used primarily as a tool for flow-field 
analysis in support of the test activities in the transonic wind tunnels of PWT. This 
paper has covered seven different problems that have utilized the PFP including both 
external and internal analysis. All but one of the examples have experimental data to 
verify the calculated flow fields, and all comparisons show excellent agreement. The 
PFP at AEDC has not been used as a tool to obtain absolute values, but rather as a tool 
to predict and verify flow fields in support of the test activities. In addition to the 
problems presented, the PFP has been used to predict the flow angularity at the model 
resulting from sting and support systems, to predict the flow around various types of 
support systems, and many other general flow analysis problems directly connected 
with wind-tunnel testing. 
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Figure 9.- Sketch of the section of the test model duplicated 
with the mathematical model. 
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Figure 14.- Upwash about cylindrical model at ^ = 0.7 and a = 2°. 
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Figure 15.- Upwash about cylindrical model with strut 
at Moo = 0.7 and a = 2°. 
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SUB-SCALE TEST DATA, AC L > Ci = 





